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ABSTRACT 

'S' : 

Q>^ . We have used the Fourier cross spectra of Cyg X-1, as obtained with BATSE 

, during a period of almost 2000 days, to estimate the phase (or time) lags between 

! X-ray flux variations in the 20-50 keV and the 50-100 keV bands as a function of 

O 

^ I Fourier frequency, v. We find that these lag spectra do not depend on source state, as 

measured by the fractional rms variations of the X-ray flux. For frequencies well below 

■ the Nyquist frequency the variations in the 50-100 keV band lag those in the 20-50 
keV band by a time interval r oc Binning effects cause the observed hard X-ray 

^ ■ time lags to decrease to zero at the Nyquist frequency. Some previous results on time 

^ : lags were affected by these binning effects. At photon energies above 10 keV the time 

', lags are approximately proportional to the logarithm of the ratio the energies of the 

^ I passbands used. 
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■ 1. Introduction 

' The delay of high-energy photons relative to low-energy photons in the X-ray brightness 

variations of Cygnus X-1 has been known for some time. The significance of these lags is their 
possible use as a diagnostic of the source of Compton upscattering of soft photons widely believed 
to be the origin of the hard spectral tail seen in Cyg X-1 and other black hole candidates (BHCs) 
(Miller 1995; Nowak & Vaughan 1996), as a probe of the dynamics of the Comptonization region 
itself (Miyamoto et al. 1988), and as a second-order statistic supplementary to the power 
spectrum in the determination of parameters in shot-noise models (Miyamoto et al. 1988). 
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The lags have been investigated using cross spectral analyses of the time variations of the 
source in different energy channels (van der Klis et al. 1995; see also Lewin, van Paradijs & 
van der Khs 1988). Miyamoto et al. (1988, 1989, 1992, 1993) have applied this technique to the 
BHCs Cyg X-1, GX 339-4, GS 2023+338, and GS 1124-338 with Ginga data, using various 
combinations of energy channels from 2-37 keV and a maximum sampling frequency of 125 Hz 
(Nyquist frequency, i^yq = 62.5 Hz). Their results showed that the lag spectra of these BHCs (i.e., 
lags as a function of Fourier frequency) were very similar, severally for the low state and the very 
high state which were covered in their observations. The amplitude of the lag in general increases 
with photon energy, E, except, perhaps, for E ^ 5 keV. This consistent picture led Miyamoto et 
al. (1992) to refer to this lag properties as 'canonical' for BHC sources. The other type of X-ray 
binary for which the lag spectra have been well studied are the Z sources. These contain accreting 
neutron stars; they show different power spectral features, and have different lag properties (van 
der Klis et al. 1987). 

In general, the results of Miyamoto et al. (1993) show a phase lag that increases as a function 
of frequency up to approximately ^^5 Hz, then decreases toward z/Nyq. At the highest frequencies, 

the uncertainties in the lag values become large enough to make the determination of any trend 
uncertain. The corresponding time lags decrease as the frequency increases towards i^Nyq where 
they are consistent with zero. At least in Cyg X-1, the phase lags depend on the average energies 
El and E2 of the energy bands used approximately as hi{E2/Ei) (Miyamoto et al. 1988). 

We have used the approximately 2000 days of data from the Burst and Transient Source 
Experiment (BATSE) on the Compton Gamma Ray Observatory to investigate X-ray lags in 
variations from Cyg X-1, in a higher energy range than previously possible. Our data cover a large 
range in hard X-ray intensity and spectral slope (see Crary et al. 1996, and the X-ray light curves 
in Paciesas et al. 1997); our observations include time intervals in which an ultra-soft component 
is present; during these intervals Cyg X-1 was likely in the high or intermediate state (see Zhang 
et al. 1997; Belloni et al. 1996). A description of the data analysis is given in the next section. In 
sections 3 and 4 we compare our results to those of Miyamoto et al., and present our conclusions. 

2. Observations 

We calculated lags between the 20-50 and 50-100 keV energy bands of the AT = 1.024 s 
time resolution data obtained with the large-area detectors (LADs). Fast Fourier transforms 
were calculated in a way identical to that discussed in Crary et al. (1996a), with the addition 
of later observations. These transforms were created for 524.288 second intervals (512 time 
bins) in two energy channels, for all detectors with an angle of less than 60° to the direction of 
Cyg X-1 (from one to four detectors view the source during a particular pointing period). The 
average number of uninterrupted 512 bin segments available with the source unocculted by the 
Earth is approximately 40 per day. The complex cross amplitudes were created from the Fourier 
amplitudes a^- = X^;;, c^. exp(z27rA;_7/n), where n is the number of time bins, c^. is the detector count 
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rate in bin k = (0 . . .n — 1) and channel number i = (1, 2), and j = {—n/2 . . . n/2) corresponds 
to Fourier frequencies 27rj /nAT. The cross spectra were normahzed to the total detector count 
rates Ni in channel i, by analogy with the method of Leahy et al. 1983), and are given by 
Cj"^ =< a'j*aj / \/NiN2 >, where the brackets refer to daily averaging. Errors on the averages 
of the real and imaginary parts of cross spectrum were calculated from their respective sample 
variances. For further averages of the cross spectrum these errors were propagated. The phase lag 
as a function of frequency is obtained from the cross spectrum as (pj = arctan [Im(C|^)/Re(C|^)]. 
With the above definitions, lags in the hard X-ray variations appear as positive angles. For these 
data, deadtime induced cross-channel effects (van der Klis et al. 1987) will affect only the real 
part of the cross spectrum, reducing it by less then 5%. 

3. Results 

3.1. Calculation of Phase and Time Lags 

Because of the properties of the Fourier amplitudes (van der Klis 1989; Leahy et al. 1983) and 
the noise in the uncollimated BATSE detectors, the cross spectra for a large number of days must 
be averaged and converted to lag values to obtain sufficiently small errors. As shown by Crary et 
al. (1996a, b) the slope of the X-ray spectrum of Cyg X-1 and the amplitude of the variability are 
strongly correlated. This indicates also that within the low state the properties of Cyg X-1 are 
determined by a simple parameter, most likely the mass accretion rate. We have therefore used 
the squared fractional rms amplitude of the noise, integrated over the 0.03 — 0.488 Hz frequency 
range, as the relevant quantity with which to correlate the lag properties. Figure Ib-d shows 
the data grouped into three parts, i.e. those for which the daily averaged squared fractional rms 
values (between 0.03 and 0.488 Hz) are greater than 0.03 and less that 0.05, between 0.05 and 
0.07, and greater than 0.07, respectively. The figure shows that to within the uncertainties, there 
is no obvious trend in the lag spectrum with source state. Figure la presents results obtained by 
averaging the cross spectra over the set of data with fractional rms squared values greater than 
0.03. 

These data show that at the lowest frequencies, the frequency dependence of the phase lag is 
consistent with zero phase lag; the phase lag rises to a peak value of 0.04 radians near 0.2 Hz and 
decreases above this point to a value to near zero radians at the i^Nyq- The resulting time lags, 
Tj = cj)j/2TTVj (where Vj is the frequency in Hz of the jth frequency bin), are shown in Figure |2| 
(filled circles). These show a roughly power-law decrease above ~0.01 Hz (exponent ~ —0.8) with 
a break above 0.1 Hz, beyond which they decrease much more quickly. 

The shape of the phase lag spectrum looks superficially like that found by Miyamoto et al. 
(1988) for Cyg X-1 and by Miyamoto et al. (1993) for GS 1124-68 and GX 339-4. For the case 
of Cyg X-1, however, the turn over in their phase lag spectrum occurred at ~30 Hz, and these 
lags decreased to near zero at their Nyquist frequency (62.5 Hz). The similarity in the shape of 
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these spectra over the samphng time dependent frequency range suggests, by analogy with the 
finite-sampUng effects in power spectral analyses (van der Klis 1989), that at least part of this 
shape (in particular the decrease toward the Nyquist frequency) is determined by binning and 
aliasing effects on the cross spectrum. 

3.2. Systematic Effects 

The effect of finite sampling on cross spectrum and lags can be calculated in a way similar to 
the analogous problem in the determination of power spectra (van der Klis 1989; see also Maejima 
et al. 1984). 

The discrete form of the cross spectrum C^^ for the pair of binned, continuous random 
functions x^{t) with zero mean (with j equals 1 or 2) in terms of the cross spectrum (7^^ of x^{t) 
and x'^{t) is given by 

fc=— oo 

where sincx = sinx/x, n is the number of bins in a data segment. At is the bin width, and 
I = —n/2 . . . n/2. The sum over k represents the aliasing in the frequency domain due to the 
uniform sampling in the time domain (van der Klis 1989), and the sinc^ terms represent the effect 
of the averaging of the data over the finite time bin width. 

To investigate the effect of binning on the phase and time lags, it is necessary to model 
the cross spectrum at frequencies above v^yq. Here we assume a cross spectrum of the form 
(712(1,) = pi^^y2nur{\u\)^ ^j^^j.^ ^.j^g modulus of the cross spectrum, is independent of the 

detailed structure of the time lag r as a function of frequency, and the difference between the real 
and imaginary parts of the cross spectrum is determined by the complex exponential depending 
on r(|z^|). This phenomenological model has the simple interpretation that for r independent of ly 
it represents a constant time lag. The real and imaginary parts of the cross spectrum are plotted 
in Figure ^. The shape of the real part of the cross spectrum is suggestive of the shape of a typical 
power spectrum, showing a 'flat top' at low frequency and a break at ~ 0.1 Hz to a power law 
form at higher frequencies (Belloni & Hasinger 1991; Crary et al. 1996b). To test the validity of 
the simple model for the cross spectrum, we have modeled F{i') with a Lor entzian- type function 
F{i') = ai/(l + (|i/|/a2)"^), and assume that the time lag has a simple power-law form, r = to|i^|~" 
from an inspection of the shape of the time lags in Figure |2| at a point well away from z^Nyq- These 
five parameters were then determined by a simultaneous fit of the real and imaginary parts of 
to the cross spectral data. The best fit parameters are found to be oi = 1.5, 02 = 0.11, 03 = 1.3 
for F(i/), and tq = 0.0071, a = -0.78 for r(|z^|). 

Using the binning correction to this cross spectrum model, the fit to the data is shown as 
the dashed line in Figure ^. Without binning, the fit form is given by the dotted line in Figure 
^. Because of the small values of r which are observed, the fit values are insensitive to the exact 
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functional from of t(|i/|), but the binning correction is a large effect. The effect of binning on the 
time lags calculated from this model is shown in Figure ^ (dashed line). We find that the time 
lag, away from i^f^yq depends on frequency, i^, approximately as up to ~ 0.2Hz. Above this 

frequency, the lags downturn, most of which is an artifact due to the data binning. 

It should be stressed that these results are not meant to provide a definitive description of the 
cross spectra, but to illustrate the important effect of binning on the lag spectrum at frequencies 
above ~ 0.5z^Nyq- In this model, the calculation of the time lags is not sensitive to the the detailed 
shape of F{v), but it is clearly not adequately described by a simple Lorentzian-like function. 
Other potential systematic effects on the cross spectrum, such as 'windowing' in the time domain 
(van der Klis 1989), have not been considered. However, the implications of this simple model 
are clear. At least part of the decrease in the time lags with frequency observed by Miyamoto et 
al. (1992, 1993) must be attributed to binning effects. To obtain a precise determination of the 
phase or time lag as a function of frequency, these effects must be considered, or, since this effect 
dominates near the z^Nyq frequency, frequencies above ~ 0.5z^Nyq should be ignored. 

Since we find that the lags are independent of rms fractional X-ray flux variations, we have 
combined our time lag results with those obtained by Miyamoto et al. (1989) for Cyg X-1. 
This allows us to study the photon energy dependence of the time lags over a much increased 
energy range. Miyamoto et al. (1989) found that the time lags between energy bands with 
average energies E2 and Ei are proportional to e = ln(£'2/i?i). For power law photon indices 
between —1.5 and —2.5 we find that for the BATSE energy bands e is in the range 0.802 and 
0.833. To avoid the binning effects we have taken the time lag at v = 0.1 Hz, for which we find 
r = 0.055 lb 0.005 s. Then the proportionality constant between time lag and the logarithm of the 
energy ratio = r/e = 0.067 it 0.006. At the same frequency we infer from Figure 3 of Miyamoto 
et al. (1989) lags of 0.045 s, 0.030 s, and 0.088 s between the 1.2-4.7 keV and 4.7-9.3 keV bands, 
between the 4.7-9.3 keV and 9.3-14.0 keV bands, and between the 4.7-9.3 keV and 15.8-24.4 
keV bands, respectively. The corresponding values of 6 are 0.045, 0.060 and 0.080, respectively 
(estimated accuracy ~ 10 percent). This indicates that between the energy ranges below 10 keV 
and above 10 keV the ratio 9 of the time lags to ln(E2/Ei) increases; however, between ~ 10 and 
100 keV 9 does not change within the accuracy with which it has been determined. 



4. Conclusions 

Based on our analysis of nearly 2000 days of BATSE data we have studied the time lags 
between the X-ray flux variations of Cyg X-1 in the 20-50 keV and 50-100 keV bands as a function 
of Fourier frequency. We flnd that the time lag spectra do not depend on source state as measured 
by the rms fractional variations of the X-ray flux. Over the range 0.01-0.2 Hz the hard lags r are 
proportional to i/"*^'*. At frequencies above O.Sz^Nyq the lags are artificially decreased (to zero at 
i^Nyq) because of data binning. At least part of the 'canonical' time lag behavior of black holes, 
as defined by Miyamoto et al., is due to this binning effect. At photon energies above 10 keV 



-6- 



the time lags are approximately proportional to the logarithm of the ratio of the energies of the 
passbands used. 
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Fig. 1. — Phase lags between 20-50 keV and 50-100 keV data for various daily averaged squared 
rms levels, s. a) All data with s > 0.03, b) 0.03s < 0.05, c) 0.05 < s < 0.07, d) 0.07 < s. 
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Fig. 2. — Time lags as a function of frequency. Solid circles denote values obtained by averaging 
over all data with fractional squared rms values greater than 0.03. The dashed line is the time lag 
obtained from a fit to the lag model including binning effects (see text) . The dotted line is a power 
law with parameters determined from the model fit. 
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Fig. 3. — Real part (upper panel) and imaginary part of the cross spectrum. The dashed line is a 
fit of the lag model to the cross spectrum with the effects of binning included. The dotted line is 
the corresponding quantity without the effects of binning considered. 



